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The infmite-C/ Anderson-Holstein impurity model is studied with a focus on the interplay between 
the strong electron correlation and the weak electron-phonon interaction. The slave boson method 
has been employed in combination with the large degeneracy expansion (1/N) technique. The charge 
and spin susceptibilities and the phonon propagator are obtained in the approximation scheme where 
the saddle point configuration and the Gaussian 1/N fluctuations are taken into account. The spin 
susceptibility is found not to be renormalized by electron-phonon interaction, while the charge 
susceptibility is renormalized. From the renormalized charge susceptibility the Kondo temperature 
is found to increase by the electron-phonon interaction. It turns out that the bosonic 1/N Gaussian 
fluctuations play a very crucial role, in particular, for the phonon propagator. 

PACS numbers: 71.27.+a,63.20.Kr,71.38.-k 



I. INTRODUCTION 

Both electron correlation and lattice dynamics are im- 
portant in understanding the physical properties of many 
condensed matter systemsi But the interplay between 
them has not been studied in detail because there are 
no reliable theoretical methods to describe the electron- 
phonon interaction in the strongly correlated electron 
systems. For weakly correlated materials we have the 
successful Migdal-Eliashberg theory^ There is, however, 
no systematic theoretical extension of the theory to the 
strongly interacting systems. With the recent advent of 
the dynamical mean field theory (DMFT) -various lattice 
problems can be mapped (or approximated) onto impu- 
rity problems. In this context it is important to under- 
stand the correlated impurity problem coupled to local 
phonons. 

One of the simplest model which incorporate both 
electron correlation and phonon is the Hubbard-Holstein 
(HH) modeli£ In the context of DMFT, it is mapped 
onto the Anderson-Holstein (AH) modelAi The AH 
model is a single-impurity Anderson model with a lin- 
ear coupling to a local phonon mode as in the Hol- 
stein models Recently, the HH model was studied with 
DMFT in combination with the numerical rcnormaliza- 
tion group (NRG) at half-filling at zero temperature^^ 
The NRG method has also been successfully applied to 
AH modelSii, and almost exact results on the electron 
and phonon spectral functions have been obtained. How- 
ever, it is still desirable to develop an analytic scheme in 
spite of its approximate nature since it helps us to under- 
stand the underlying physics in more clear and intuitive 
way. In this paper, we will consider the regime where the 
electron correlation dominates over the phonon effect, 
which will allow us to study the phonon effects on the 
Kondo physics. The pure Holstein model without elec- 
tron correlation is studied in Ref. ^3 using semiclassical 



approximation and in Ref. [Tl| using NRG and DMFT. 

When the local Coulomb repulsion of Anderson model 
is very large compared to the hybridization between the 
impurity and the conduction electrons, the impurity de- 
velops a magnetic moment at high temperature regime^ 
The moment is gradually screened by the spin of the con- 
duction electrons as the temperature is lowered below 
the Kondo temperature scale. The physical properties of 
Anderson model at temperature lower than the Kondo 
temperature can not be described within perturbative 
framework since the effective Kondo coupling constant 
becomes very large in that temperature regime. This low 
temperature state is often called the Kondo fixed point 
in the sense of renormalization groupAi^ We mention 
that the physics of Kondo fixed point can describe the 
correlated metallic states such as heavy fermions in the 
context of DMFT. 

In the limit of the infinite local Coulomb repulsion the 
impurity develops the local moment at high tempera- 
ture and the system crosses over into the (completely 
screened) Kondo fixed point at low temperature. Thus 
the physical properties of the system in the limit of in- 
finite local Coulomb repulsion are expected to be qual- 
itatively similar to those of other physical systems with 
very large but finite local Coulomb repulsion. If the lo- 
cal Coulomb repulsion is taken to be infinite the dou- 
bly occupied state of the impurity is not allowed in the 
physical Hilbert space. In general, the standard per- 
turbative approach based on the Wick theorem on the 
electron operators and the associated Feynman diagram 
expansions does not apply to the systems with the con- 
strained Hilbert space. This fact poses a technical diffi- 
culty for the approaches based on the electron operators. 
However, if we employ the slave field representation of 
electron operators, the constrained Hilbert space can be 
treated in a considerably simplified way. The slave field 
method also admits the Feynman diagram expansion and 



2 



the direct manipulations of functional integrals. This 
slave field method has been successfully applied to the 
infinite-t/ Anderson modeliii^i^ and it is particularly 
adequate in the study of the low temperature properties 
of the Anderson model. The bose condensed state which 
emerges in the saddle point (mean-field) approximation 
embodies the Kondo resonant state, which is the crux 
of the Kondo fixed point. Thus the essential physics of 
the Kondo fixed point is captured in the leading approx- 
imation of the slave boson method. At high temperature 
the bose condensate disappears and it corresponds to the 
high temperature weak Kondo coupling regime where the 
perturbative approach such as the poor man's scaling 16 
can be applied. 

Our objective in this paper is to incorporate the 
phonon dynamics into the infinite- U Anderson model 
within the scheme of the slave boson method along the 
line pursued by Coleman^ We will assume that the 
electron-phonon interaction is weak compared to the 
strong electron correlation, so that the physical situa- 
tion at the zeroth approximation is still best described in 
terms of Kondo fixed point physics. We have developed a 
scheme in which the effect of electron-phonon interaction 
can be computed systematically in the large-degeneracy 
limit or, in other words, the l/N expansion. The basic 
idea of the l/N expansion is briefly illustrated in Ap- 
pendix E] 

The charge and spin susceptibilities of the impurity 
and the local phonon propagator have been obtained. 
Note that the dynamical correlation functions are gener- 
ally very hard to compute when the electron correlation 
is very strong. The feasibility of the analytic computa- 
tion of the dynamical correlation functions is the great 
advantage of the slave boson methods over other meth- 
ods. It is found that the spin susceptibility of the im- 
purity [Eq. (|67fl ] is not modified by the electron-phonon 
interaction up to the approximation of Gaussian fluctu- 
ations. The charge susceptibility of the impurity up to 
the second order in the (weak) electron-phonon coupling 
is given by Eq. (|77J) . From the charge susceptibility at 
zero frequency, the Kondo temperature renormalized by 
the electron-phonon interaction can be extracted, and 
it is given by Eq. Q81[l. Finally the phonon propaga- 
tor is found to be Eq. 1)82(1 . These are the main re- 
sults of this paper. In all of these results the influence 
of electron-phonon interaction is small since the charge 
fluctuations which directly couple to the local phonon 
are severely suppressed by very strong electron correla- 
tions. These results manifestly demonstrate the essen- 
tial interplay between the strong electron correlation and 
the (weak) electron-phonon interaction. For instance, the 
small renormalization by electron-phonon interaction in 
Eq. (|81|l is expressed by a product of the polaron en- 
ergy and the Kondo energy, where the nonpeturbative 
physics of the electron correlation is encapsulated in the 
Kondo energy scale. The smallness of the influence of the 
electron-phonon interaction in our case is partially due to 
the absence of the orbital degrees of freedom. The effect 



of phonons can be much more substantial for the materi- 
als with multi-orbitals, where the phonons can couple to 
the orbital quantum number^ Our approach can be ex- 
tended to such cases, and the progresses in this direction 
are being made. 

This paper is organized as follows. In Sec. II we intro- 
duce the AH model and the slave boson formalism. Then 
the rescaled effective action in the Read-Newns gauge is 
obtained. In Sec. Ill we summarize the basic results on 
the saddle point approximation by Colemanpii In Sec. 
IV the phonon dynamics is coupled to the bose excita- 
tion emerging in the 1 /N expansion and the effective ac- 
tion of them is computed. In Sec. V the spin and the 
charge susceptibilities are calculated from the effective 
action obtained in Sec. IV and their physical properties 
are analyzed. In Sec. VI the phonon propagator is cal- 
culated from the effective action obtained in Sec. IV. 
Sec. VII is for the summary and concluding remarks. 
Some of details of the calculations can be found in the 
Appendices. 

II. FORMALISM 
A. Slave boson representation 

The Hamiltonian of AH model consists of three parts. 

UAH =T~iel +T~iph+'Hel-ph- (1) 

H e i and 7i p h is the Hamiltonian for electron and phonon 
part, respectively, and H-ei-ph is the Hamiltonian describ- 
ing the interaction between the electron and the phonon. 




H et - ph = g Q(£ l tif<r -(£,&*))■ ( 4 ) 

a cr 

a =t, I is the spin index, and Vk is the hybridization ma- 
trix element. f a is the impurity electron operator, and 
Cka is the conduction electron operator, e/ is the energy 
of the impurity level, and U is the local Coulomb repul- 
sion at the impurity site. Q is local phonon coordinate. 
P is the conjugate momentum of Q satisfying [Q, P] = ih. 
M is ion mass and fl is the oscillator frequency of disper- 
sionless (Einstein) phonons. iV) at is the number of lattice 
sites for the conduction electrons, go is the (unsealed) 
electron-phonon coupling constant. 
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In the limit of infinite U, the doubly occupied state 
of impurity / electrons is prohibited. In this constrained 
Hilbert space, the impurity electron operator f a can be 
expressed in terms of fermion operator s a and the slave 
boson operator b satisfyin g 13 ' 17 



(5) 
(6) 



fa = 4 b ) U = s <? bt > 



Now the remaining empty state and the singly occupied 
states denoted by spin quantum number a at the impu- 
rity site can be represented as follows. 



| empty) = 6 T |ref), 
\a) = 4|ref), 



(7) 



where |ref) denotes some arbitrary reference state. The 
local Coulomb repulsion term U /j j "j/j which is non- 
zero only for the doubly occupied states can be dropped 
in the constrained Hilbert space. After the elimination 
of the doubly occupied state we can formally increase the 
number of spin components from 2 to arbitrary N. 

The Hamiltonian in the constrained Hilbert space can 
be expressed in terms of the slave fields as follows, (the 
index a has been changed to m in anticipation of l/N 
expansion) 



N 



k,m—l 

i / i 



lai 



P< 



-(mq 2 q 2 . 

2V M > 

9oQ(%2 4l S m - 



(8) 



The electron-phonon coupling gives rise to the charge 
fluctuation since it couples to the total charge at the 
impurity site which is cither or 1 in the constrained 
Hilbert space. Note the scaling of the hybridization ma- 
trix element V k — > \J^V k in Eq. (JSJ. /i c is the chemical 
potential for the conduction electrons. At this point we 



pass into the Lagrangian formulation in imaginary time. 
L = >^ c\ m d T + e k - fi c c km 

km 

+ Y 4 [ d r + e /] s m + tfd T b 

m 

■((a r Q) 2 + n 2 Q 2 



VMat 

M 

Y 

9oQ(J2 



\C£sls m + tfb-l), 



(9) 



where A is a Lagrange multiplier implementing the con- 
straint Eq. 10. In the l/N expansion this constraint 
will be realized through the conserving property of the 
l/N expansion^ The partition function is given by the 
following functional integral 



Z = J D[b, A] J D[s,Q] J D[c]e 



f cItL 



(10) 



The integration over the conduction electrons c km can 
be done exactly and it generates a self-energy Eo(ie) for 
the slave fields. 



E (ie) 



— V- 



Niat ^ ie + [i- e k 
ReS - »A sgn(e), 



A = nN(E F )(\V k \* 



>FS, 



(11) 



where N(Ep) is the density of states of the conduction 
electrons at Fermi energy and the subscript FS denotes 
an average over Fermi surface. The hybridization matrix 
element V k was assumed to be weakly dependent on mo- 
mentum near the Fermi surface of conduction electrons 
in the second line of Eq. (fill) . From now on ReEo will be 
neglected since it is the order of the inverse of energy cut- 
off. After the integration over the conduction electrons 
Cfem, the resulting action becomes 



S = 



/ dr V 4n( 5 r + £/)s m + 6 f 5 T 6 

+ 1 [\rdr'J2^o(r-r')sl(T)b(r)s m {r')b^T') 

+ jyrf[(d T Q) 2 + n 2 Q 2 ) 
drg Q 

m m 

/ drA(V;4 Sm + 6t6-l), 
Jo 



(12) 
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where 



S (r-r') = T^ e --(— ')s (ie). (13) 



B. Read-Newns gauge and Rescalings 

The systematic computations of 1/N fluctuations can 
be most clearly achieved in the Read-Newns (RN) 
gauge 



13.14.15 



b(r) = R(t) e l6 ^\ s m (r) = z m {r) 



(14) 



Note that the original /-electron operator is invariant 
under the following gauge transformation: 



s m e lv , b — > be 



(15) 



R(t) > is the modulus of complex boson b. Following 
Coleman^ we write the fluctuating Lagrange multiplier 
as 



A = if2 + Xsa, 



(16) 



where X sa is the real part of the saddle point value of 
the Lagrange multiplier, which is to be determined later. 
In RN gauge, the infrared divergences which plague the 
expansion in terms of the original fields 6, s m can be 
avoided. 

The action in RN gauge is expressed as follows: 
rP 

S= d,T^zj n (d T + e f )z m 

+ [ drdT'J2Mr-T')zl(T)R(r)z m (r')R(r') 

+ f dT Y( {dTQ)2 + n2Q2 ) 

+ / drgoQ 

z m z m . z }n z m) ) 

f dO 

+ dTi{—+VL )C^z] n z m + R 2 -qN) 

•'° T m 

[P 

+ / dr X sa (R 2 - qN), 
Jo 



where 



Se f = e f + X sai q=l/N. 



(17) 



(18) 



The integration of phonon Q contributes the following 
term to the action 

2 1-0 rP -iw(t-t') 

/ ^ / dr'TT 6 . n2 
x^<5n sm (r)^5n sn (r'), (19) 



where 5n sm = sjj s m — { s m s m}- in order for the 1/N 
expansion to be applicable, each term of the action should 
be proportional to N. The structure of Eq. (fTTif suggests 
that the electron-phonon coupling go should be rescaled 
as follows: 



go = g/VN. (20) 
By the same reason R(r) should be also rescaled. 

R( T ) = VNr{T). (21) 
The rescaled action becomes 

rP 

+ / dTdT'J2Mr-T')zl(T)r(T)z m (T')r(T') 

+ f dr^({d T Q) 2 + tfQ 2 ) 

If 13 /' 
V m m 

+ f dTi(^+n )^zlz„ 
r 13 do 

+N / d Ti{ — + n )(r 2 -q) 



rP 

+N / dT\ sa (r 2 -q), 
Jo 



(22) 



Now the 1/N expansion can be directly applied to the 
action Eq. 



III. SADDLE POINT APPROXIMATION 

The leading order approximation of the 1/N expan- 
sion amounts to the saddle point approximation or equiv- 
alently, the mean-field approximation. The mean-field 
theory of infmite-i7 Anderson model has been studied in 
great detail in Sec. II of Coleman's paper^ Here, we 
briefly summarize the main results of Ref^ for the pur- 
pose of the self-contained treatment and the introduction 
of necessary notations. 

We assume that r, A take a static value r sa , X sa in the 
saddle point approximation, respectively. In the saddle 
point approximation the following flucutations with re- 
spect to the saddle point value 



d6 

i(— + Clo), 5r = r-r s 



(23) 



are neglected. We will treat the electron-phonon inter- 
action on an equal footing with the above fluctuations, 
so that the electron-phonon interaction is not included 
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in the saddle point approximation. Now the action ap- 
propriate for the saddle point approximation can be ex- 
pressed as 

rP 

S m f = dT Y] z L( d r + e f )z m 

+ [ drdr'VEo(T - T')r 2 8a zl(r)z m (r') 
Jo 



+N / drX sa (r 2 a -q), 



(24) 



Next we integrate over z rn to obtain the saddle point 
action S sa (r sa , X sa ) which is 



(25) 



q 


A/A 


e//A 


A A -/Ao 


1/2 


2.80 x ltr- 5 


2.46 x l(T b 


2.80 x 10~ d 


1/4 


1.95 x l(T a 


1.99 x lCT 3 


2.79 x 10~ 3 


1/6 


1.39 x 10" a 


2.41 x lO" 3 


2.78 x 10- a 



TABLE I: The numerical solutions of the saddle point equa- 
tions Eqs. 1)291310 . The input parameters are D = 65 Ao and 
t f = -3.2A . A K is defined in Eq. <|3H . 



where 



•*A, 



A 



In 



nD 



(33) 



The explicit form of S sa is given by 

S sa = N dr X S a(r 2 sa - q) 
Jo 

- iVTrln[a r + ^+r s 2 a So]. 



(26) 



The saddle-point values r sa , A sa are determined by mini- 
mizing S sa with respect to them. The minimization with 
respect to A sa yields 



- T E — 



1 



e - e/ - A sa + iAsgn(e) ' 



where 



A = A r s 2 Q . 



(27) 



(28) 



At zero temperature the summation of Eq. (|27|1 can be 
done in a closed form. 



q 



2 7T V A ' 



(29) 



Likewise, the minimization with respect to r sa yields 

-iA sgn(e) 



A., 



ie - e/ - A sa + iAsgn(e) ' 



(30) 



At zero temperature the summation of Eq. (|30|l can be 
done in a closed form. 



X, 



Ac 
2vr 



In 



D 2 



Ai + (e f + X sa y 



(31) 



where D is an energy cutoff which is the order of the 
bandwidth of conduction electrons. Two real-valued 
equations Eq. Q29H and Eq. H31|l can be combined into 
a single complex-valued equation. 



«A g, 



(32) 



The Kondo limit is specified by the condition |£| <§C Ao, 
and in this limit the first term of the right hand side 
of Eq. (|32|) dominates. In the Kondo limit, the Kondo 
temperature scale is given by 



lei 



e 2 ~De^ /Ao . 



(34) 



This result derived from the saddle point approximation 
coincides with the one obtaind from the scaling theory of 
Anderson model by Haldane except for the prefactorii§ii£ 
A few numerical solutions of Eqs. (|29l31|l are presented in 
Table I. From the numerical solutions we find that in the 
case of q = 1/2, A 3> e/, while for q = 1/4 and 1/6, 
A ~ Ef. In all cases Ak and A are of the same order of 
magnitude. 



IV. 1/N FLUCTUATIONS 

The 1/N corrections are related to the fluctuations 
with respect to the saddle point configuration. It is con- 
venient to treat the fluctuations collectively so let us 
define a two-component vector X as follows: 



X 

e 



dr 



d9 

— +n . 

dr 



(35) 



The 1 /N corrections can be computed systematically by 
expanding the action Eq. (|22|) with respect to the fluc- 
tuations and by integrating out z m . The action Eq. I|22|) 
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can be re-expressed as follows. 

S = S% + S z + S p h + SS Z , 



s z 

Sph 

6S Z 



^ [ — ie + if — iAsgn(e) 

SI 



In 

10 v 



(36) 
(37) 
(38) 

(39) 



dl~ r-^{ z ln z rn) 

V JV 
r P 



rP rP 

x (r sa £r(r') +r sa 5r{r) + (5r(r)fr(r')) , (40) 



where 



r<°> = 2 



Asa 1 

1 o/ ■ 



(41) 



<f> m is the source field necessary for the computation of 
impurity susceptibilities. The irrelevant constants are 
omitted in Eq. 

Next step is to integrate over z m . This step involves 
S z and 5S Z of the action Eq. i|36|) . To simplify notations, 
define 



G 1 (ie) = ie + iAsgn(e) — if. 
The integration over z m generates 

- ^ Tr In [ - G (0) (r - r') + M (r, r') 

m 

where M(r, r') is given by 



M(r, r') = (0 m + -£=Q + »0) (r)i(r - r') 

v iV 



(42) 



(43) 



+E (r - r')(r sa Sr(T) + r sa 5r{T r ) + 5r(r)5r(r')j(44) 

In our approximation we expand the logarithm of 
Eq. (|43|l with respect to M(r, r') up to the second order. 
This is equivalent to reckoning in the Gaussian fluctua- 
tions but not the higher order fluctuations in 1 / >/~N (see 
Appendix \K$ . 

The relevant part of the expansion is 

+^tt[g m] +^J2[g mg m . 

m m 

The details of further calculations are presented in the 



appendix The effective action of X, Q, 4> m is given by 



N 



, /7 T E^Nr (1) ^H 

iuj 

x (a^ } (z0r SQ + <£ m r so ) + *r) (-iw) 
+ ^ee'C^) {i0r sa ){iu) (0 m r sa )(-iw) 

m iw 

sa 0m 

(— ia;) 

m ia; 



where 

r (1) (iw) 



2 A 4- V /T W 2 + ft'' 



2 + K. 



(i) 



A'. 



(i) 



(46) 



The explicit forms of various polarization functions K 
can be found in Appendix [U] By the saddle point condi- 
tion Eq. lHU, 



2A sa + K<$ = 0. 



Then it follows that 

r (1) M = 



2 + K, 



(i) 



K, 



(i) 
oo 



(47) 



(48) 



In fact, only Re K$ in the (1,1) element of Eq. 
contributes to the action due to the symmetry lu <-> — lj. 
Now the integration over phonon 

exp[-S$ f ] = J D[Q] eMSiff(X,Q,c(> m )] (49) 
contributes the following to S^ff- 



E, 



P 2 
T 



n 2 



(K$r sa (iG + <f> m )(ico) + K$8r(w) 



where 



(,50) 



(51) 



is the polaron energy. Since we are assuming that the 
electron-phonon correlation is weak, the relation 



E p < A 



(52) 
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holds. We will also assume that E p > Ak, which is 
more relavant to the real physical systems. The effect of 
electron-phonon interaction can be traced by following 
E p . Let us define the following dimcnsionless function 



D(iu) 



n 2 



+ nz + K$g*r* a /M 



(53) 



The explicit form of the effective action S^^{X, <j> m ) is 
given by 



m.n.iuj 



(1 - E p D(iu;y sa K { e 

m iuj 



where 



T{iuj) = 
K rr 



K r 6 Keg 

Kg - E p rl a D{iu)[K$]\ 
2 + K$-E p rl a D(iu > )K$Kg\ 
^(l-^iw)^). (55) 
Finally the integration over X 

expI-S^ = J D[X] exp[-Si%(X,4> m )} (56) 
yields the effective action of <j) m . 

S e C ff = \^^ K w( iu} ) r2 sa ( t ) rn{iu)(t> m {-iuj) 
m ioj 

E p r 4 sa ^ [K^Dii^^LO^i-tLu) 



2N 



m.n.iuj 



(ih 2 



X ^(l^M-l^lKi^K^T- 1 



(57) 



V. THE SPIN AND THE CHARGE 
SUSCEPTIBILITIES OF IMPURITY 

The spin susceptibility of impurity is defined by 

Xs(t-t') = (M z (t)M z (t')), (58) 
M z = 2J ms] n s m = ^ mz m z m, (59) 



where the sum over m runs from —N/2 to N/2. The 
charge susceptibility is defined by 



where 



Xo(r-T') = (5n f (r)5n f (r')), (60) 



5fi s = E s ™ Sm ~ (E s ™ Sm ) 

m m 

= E z ™ z ™ ~ (E z ™ z ™) ( 61 ) 



The spin and charge susceptibilites can be computed 
from 



Xmn ( 7_ ; T~ ) 



S 2 InZ 



50 m (r) <5</>n (r') 



6—0 



(62) 



Here the partition function Z is understood to be a func- 
tional of 4> m only, namely all other fields should be inte- 
grated out. In our approximation, 



Z = J D\<j> m ] exphSg 3 /^)], 



(63) 



where S^) f {4> m ) is given by Eq. Xmnfar') in frc- 



quency space is given by 

Xmn(r,T') =TJ2 Xmnii^e-^-^. 

Then it follows that 

Xs(iu) = mnxmn(iuj). 

m.n 

From Eq. we find 



(64) 



(65) 



Xr. 



r -^{l-E p D(iu>y sa K$) 2 



-fE p D^) [ffgW* 



(66) 



From Eq. I|66[l the spin susceptibility is readily obtained 
Xs(iw) = [E™ 2 ] X s ( iuj ), 

m 

X.(iw) = -r 2 sa K$(iu). (67) 

In our approximation scheme the spin susceptibility is not 
renormalized by the electron-phonon interaction. Clearly 
this is an artifact of the approximation of the saddle point 
and the Gaussian fluctuation. If the fluctuations be- 
yond the Gaussian fluctuations ( higher order in l/y/N) 
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are considered the spin susceptibility is expected to be 
renormalized by the electron-phonon. However, the cor- 
rections due to higher order fluctuations will be rather 
small. It is illuminating to compare this result with that 
of numerical renormalization group studies. Hewson and 
Meyer have computed the imaginary part of the retarded 
spin susceptibility for the symmetric and finite-U AH 
model (see the Fig. 13 of Ref. |Q. Compare two cases 
with A = 0.0 and A = 0.02). Their result shows that the 
electron-phonon interaction does not influence the qual- 
itative behavior of the spin susceptibility except for the 
decrease of the peak height at very low energy around 
the Kondo temperature scale. For our parameter regime 
where the electron correlation is far more dominant over 
the electron-phonon interaction, the decrease of the peak 
height would be hardly noticeable. 

The charge susceptibility X c{iw) is given by 

Xc(ioj)/N = r 2 sa (l - E p D{iu)r 2 sa K$) 2 



-r z K 

' so 



${iw)(l-E p D{iu)r 2 sa K { 



Xc(iu). 



(68) 



The charge susceptibility is indeed renormalized by the 
electron-phonon interaction in a very complicated way. 
This complicated expression can be recast in a more 
transparent form in the following way. First note that 
the charge susceptibility in the absence of the electron- 
phonon interaction is given by 



(1) Jrii'l 



[r(i) 



(i) 

r0 
(1) 



K, 



(i) 



(69) 



This result can be obtaind from Eq. (|68|l by putting E p = 
everywhere and was derived by Coleman. 13 Explicit 
calculation shows that 



K, 



(i) 



(^ } +2) 2 



(70) 



For the subsequent analyses it is helpful to examine the 
susceptibilities in the absence of phonons in more detail. 
The susceptibilities at zero frequency in the absence of 
phonons are given by 



xi o) (^ = o) = --=-, 



1 A 



xi o) (^ = o) = 



ttA 



(A + 7r? / ) 2 + ( 7 rA) 5 



(71) 



In the Kondo regime, the relation Ao tt^/A 
holds, so that we have 



TTA 



A 2 



~2 



(72) 
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FIG. 1: The real and imaginary part of the retarded charge 
susceptibility in real frequency in the absence of electron- 
phonon interaction (iu) — > uj + iS). The dashed line 
and the solid line indicates the real part and the imaginary 
part, respectively. Note that the frequency is in the logarith- 
mic scale, and the y-axis is in the absolute scale. The input 
parameters are identical with those of Table I. The value of q 
is 1/2. 



The ratio of the charge and spin susceptibility at zero 
frequency is then given by 



X (°>^ = 0)/xi 0) (^ = 0) 



ttA 



K 



An 



< 1. 



(73) 



In fact, over most of the frequency range of interest, the 
absolute magnitude of spin susceptibility is far greater 
than that of the charge susceptibility. This is because the 
charge fluctuations are severely suppressed by the very 
strong local Coulomb repulsion in the low temperature 
Kondo regime, while the spin fluctuations (spin flips) are 
the dominant processes there. 

We also note that is related to the propaga- 

tor (5r(ioj)5r(—iuj))( ) in the absence of phonon in the 
following way 



(frHM- M ))(») = ll 



(74) 



which can be shown by inverting the matrix T^iiio). 
This implies that higher order corrections stemming from 
the bose fluctuations are suppresed in the limit of large 
N. 

Using the above results on the susceptibilities in the 
absence of phonon and expanding the matrix up to 
the second order of g, the result Eq. I|68l) in the presence 
of phonon can be expressed in a more transparent form: 



U^) * X ( c ] M [l + E p D{ilu) xf ] M] , (75) 

where D(iuj) is defined by Eq. (|53() . Using the previous 
results we find 



EpD^x^iiLu) < 1. 



(76) 



9 
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FIG. 2: The imaginary parts of the retarded charge suscepti- 
bility and spin susceptibility in the absence of electron-phonon 
interaction. The charge susceptibility has been multiplied by 
a factor of 25,000 to make comparison easier. The input pa- 
rameters are identical with those of Table I. q = 1/2. 



FIG. 3: The real parts of the retarded charge susceptibility 
and spin susceptibility in the absence of electron-phonon in- 
teraction. The charge susceptibility has been multiplied by 
a factor of 10,000 to make comparison easier. The input pa- 
rameters are identical with those of Table I. q = 1/2. 



The structure of Eq. (|75(l clearly prompts us to write 



~(0)/. \ 
Xc 



(77) 



The imaginary part of the retarded charge susceptibility 
[Eq. (|75J) ] in real frequency can be expressed as 

ImxfH ~ Imx^^'M 

x (l + E p Re[D R (w)}Re[x^ 0) (")}), (™) 



where only the leading contribution is exhibited. As de- 
picted in Fig. ^ the variations of Re[xf ^ (u)] is almost 
negligible for all frequency region below the frequency 
of the order of |e/|. Therefore, the major frequency de- 
pendence of the electron-phonon correction to the charge 
susceptibility originates from Re[D R (oj)]: 



Re[D R (w) 



n 2 



ft 2 ' 



(79) 



where the singularity at w = ft is cut-off by the imaginary 
part of Kqq but this is irrelevant in our discussion. In 
view of Eq. (|79|l a qualitative change near the phonon fre- 
quency ft is expected. More specificially, one can expect 
that the imaginary part of the charge susceptibility is to 
be enhanced at frequency slightly lower than the phonon 
frequency while it is expected to be reduced at frequency 
slightly higher than the phonon frequency. This expecta- 
tion is in good agreement with the results of NRG studies 
by Hewson and Meyer (see the left panel of Fig. 13 in 
Ref. and compare two cases of A = 0.0 and A = 0.02). 

The charge susceptibility at zero frequency in the pres- 
ence of phonon can be expressed 



7rA 



ttE„A 
A 2 



(80) 



For the case of q — 1/2 the relation A^- ~ A holds, and 
the result Eq. I|80|) can be interpreted from the viewpoint 



of the renormalization of Kondo temperature by electron- 
phonon interaction. 



T, 



K 



T, 



(o) 



K 



A 2 



il 



where T% — Ak is the Kondo temperature in the 
absence of phonons. The Eq. I|81J) indicates that the 
electron-phonon interaction increases the Kondo temper- 
ature, which is consistent with the recent results found 
by NRG methodM. Since E p < A and A < A the 
amount of increase is rather small. This smallness of 
renormalization is again due to the suppressed charge 
fluctuations. With Coulomb repulsion fixed to be infi- 
nite from the outset, the charge fluctuation is controlled 
by the e/. As the impurity level goes down deeper into 
the Fermi sea of conduction electron, evidently the charge 
fluctuation is reduced. As can be seen in Eq. I|34(l . A de- 
creases exponentially with the decreasing e f (in the sense 
of becoming more negative) . 



VI. PHONON PROPAGATOR 

The phonon propagator can be obtained from the ac- 
tion Eq. by setting the source field <fi m to zero and 
then by integrating over the bose fluctuations X and the 
fcrmions z m . 



D ph (ioj) = \mlo 2 + Afft 2 - 5 2 xi 0) M 



(82) 



The physics of Kondo fixed point is built in Xc\i^)- If 
the bose fluctuations X were not taken into account, the 
phonon propagator would be 



Mlo 2 + Mft 2 - g 2 xi 0) (io;) 



(83) 



This can be seen from Eq. I|45|l by setting the bose fluc- 
tuation X and the source field <p m to zero, and then by 
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integrating over z m . As discussed in Sec. Eland can be 
seen Fig. |2]and Fig. [3]tlie charge and spin susceptibility 
of impurity behave essentially differently from each other. 
First of all, the absolute magnitude of the spin suscepti- 
bility is much larger than that of the charge susceptibility 
by the order of magnitudes. Second the frequency depen- 
dences are also different from each other. The imaginary 

part of the spin susceptibility Imxf''°'(w) has a sharp 
peak at to ~ , while the imaginary part of the charge 

susceptibility Im^f ' (u>) has a rather broad peak at fre- 
quency higher than A^. Imxf''°'(w) also has a peak at 
high frequency order of — e/ as seen in Fig. ^ 

If the bose fluctuations were not considered, then based 
on Eq. H83f) one would anticipate rather substantial renor- 
malization even for the relatively small values of g. This 
is in stark contrast with the results by NRG 6,7 which 
show rather small renormalization. To make arguments 
more concrete let us compare last two terms in the 
bracket of Eq. (|83[) at zero frequency. Using the result 
Eq. JTlJ we obtain 



2~(0)/- 

9 Xs («w 



0) 



Mfl 2 



3l >1 

A k 



(84) 



for physically relevant cases. 

However, regarding the result Eq. i|82|l for which the 
bose fluctuations taken into account, the ratio becomes 



2~(0)/- 
1 Xc («w 



0) 



mq 2 



E p Ak 

A2 
^0 



«1, 



(85) 



which agrees with the results by NRG. Thus, unless the 
electron-phonon interaction is extremely strong, the self- 
energy correction to the phonon propagator coming from 
the electronic sector is small. But if the electron-phonon 
interaction becomes so strong that the ratio Eq. I|85l) be- 
comes order of unity, the saddle point we have started 
with won't be valid any more and new saddle points will 
emerge. The results by NRG studies indeed demonstrate 
that new features emerge in the strong coupling regime of 
electron-phonon interaction. Unfortunately our approx- 
imation scheme does not apply in that regime. In any 
case we conclude that to obtain the physically correct 
result it is absolutely essential to include the 1/N bose 
fluctuations X. 

The form of the phonon propagator Eq. I|82l) suggests 
that the phonon frequency softens. The renormalized 
phonon frequency is given by 



_ 7T EpA 

2 A 2 



(86) 



The result Eq. I|86|l can be compared with that de- 
rived in the weak (electron-electron and electron-phonon) 
coupling random phase approximation type perturbative 
calculations^ 



V TrAn/ 



(87) 



The result Eq. 1|86[) should be also compared with that 
derived in the semiclassical approximation for the large 
but finite-U and symmetric AH models 



n 1 



EpAo 
~U 2 ~ 



The results Eqs. Ij86l88(l indicate that the strong Coulomb 
repulsion suppresses the softening of phonon frequency 
but the detailed form of the softening depends crucially 
on the details of Hamiltonian such as the symmetric or 
asymmetric nature of charge fluctuations of Anderson 
model. 



VII. SUMMARY AND CONCLUDING 
REMARKS 

We have studied the AH model with the infinite lo- 
cal Coulomb repulsion. The limit of the infinite local 
Coulomb repulsion eliminates the doubly occupied states 
in the Hilbert space of impurity. The slave boson method 
is very effective in treating the problems with the con- 
strained Hilbert space, where the strong correlation is ac- 
tually built in the constraint itself. In general the slave 
boson method is more reliable in the low temperature 
regime where the thermal fluctuations are small. The 
1/N expansion scheme is employed to control approxi- 
mations systematically and it is found that the electron- 
phonon interaction can be naturally fit into the scheme. 
Remarkably the nonperturbative physics of the low tem- 
perature Kondo fixed point is captured in the leading 
approximation. We have used the saddle point approx- 
imation supplemented with the Gaussian fluctuations. 
In fact, the Gaussian fluctuations play a very crucial 
role in determining the physical properties of impurity 
and phonon. The spin susceptibility is not renormal- 
ized up to this approximation. The Kondo temperature 
is found to increase by the elctron-phonon interaction 
but by a rather small amount. The same is true of the 
charge susceptibility and phonon propagator. This is be- 
cause the quantum fluctuations of phonon which couple 
to the charge fluctuations of the impurity electron are 
suppressed by the very strong local electron correlation. 
But the interplay between the strong electron correlation 
and the electron-phonon interaction is clearly displayed. 

We could also have calculated the electron spectral 
function but all the important features such as the in- 
crease of the Kondo temperatures and the pronounced 
phonon effect near the phonon freqency are essentially 
identical with those of the charge susceptibility. 

The AH model can be naturally generalized to the im- 
purity models with the orbital degrees of freedom. As 
a matter of fact, this is much more realistic model than 
the AH model since the electron correlation effects are 
more prominent in the d and / electron systems where 
the orbital degrees of freedom are active. For these cases 
the simple Holstein phonon should be also generalized to 
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incorporate Jahn- Teller type phononsi Slave boson ap- 
proach can be extended to the orbitally active cases, and 
currently the progresses in this direction are being made 
by the authors. 
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of Ref. 0, and here we give an equivalent but techni- 
cally slightly different treatment. The measure of boson 
functional integral is given by 



D[r,0]=H 



d(r 2 (r)9(r)) 



2tt 



(Bl) 



The fermion measure is invariant under the RN gauge 
transformation. 



D[a ) a*] = D[z,z i ]. 



(B2) 



The angular variable 9 in the sector with the winding 
number m can be expanded 13 



)(r) = 



+ 27rm-, (B3) 

P 



APPENDIX A: ILLUSTRATION OF THE BASIC 
IDEA OF 1/N EXPANSION 

Let us briefly illustrate the basic idea of 1 /N expansion 
in the case of one-dimensional integral. Consider an inte- 
gral I of the following form defined on a certain domain 
D: 



I = I e~ Na{x \ iV> 1. 

ID 



(Al) 



If a minimum of g(x) is attained at xq € D, then the 
function can be expanded 



g(x) ~ g(x ) + l-g {2) (x )(x-x Q ) 2 + 



3 (n) (^o) 



n>3 



(x-x ) n , 
(A2) 

where denotes the ra-th derivative and g( 2 \xo) > 0. 
Defining a rescaled variable z 



(A3) 



the integral I can be approximated by 

I ~ e -Ng(x )-lnlNg (2 > (x )]/2 / dz £ -z 2 /2- Cl z 3 /VN+- 



^N g (x )-\n[N g ^(x )]/2 ( , COnSt - , . \ M) 



const. 
N 

The integration over z can be done systematically as a 
power series expansion in 1 / \fN . In this paper we ignore 
the corrections of order 1 / y/N and higher in the bracket 
of Eq. (|A4|) , but they can be computed in a controllable 
way if desired. 



APPENDIX B: DETAILS OF THE 
CALCULATIONS OF 1/N CORRECTIONS 

First let us specify the measure for the functional in- 
tegral in the RN gauge. This is discussed in Sec. Ill B 



where v n is the bosonic Matsubara frequency given by 
v n = 2im/p. This expansion automatically satisfies the 
winding number condition without further conditions im- 
posed on the coefficients 9 n . The zero mode term 9q 
should be factored out since our action does not depend 
on it. For the variation which involves only the zero 



mode we have d9^' l> — l/\/j3d9o. Factoring out J d9 { 



■) 



we obtain 



Dm = e ^ n 



d9n 

2tt 



(B4) 



Once the zero mode term is factored out both 9(t) and 
d9(r)/dT are characterized by the same set of coefficients 
{#n^o} and an integer m. Explicitly 



d0 n 

— = 27TTO (3- 

dr 



From the above expression we find that the constant part 
of is constrained to take a value which is an integral 
multiple of 2tt/(3. However, as can be seen in Eq. ill' 21) 
4^ always appears in combination with f2o- Since f2o is a 
part of Lagrange multiplier implementing delta function 
constraint it is non-compact, in other words, it can take 
an arbitrary real value. Now the explicit summation over 
the winding number is unnecessary and we need only to 
consider a bose field 



9 = | + tt , 



(B5) 



which has an ordinary mode expasion without any con- 
straint. Effectively we can replace D [9] with -D[0]. 

If the fluctuations beyond the Gaussian fluctuation of 
the order 1/N can be ignored, namely the corrections of 
order 1/N 3 / 2 and higher can be ignored, the functional 
measure for the r integration can be approximated by 



dr 2 {r) = 2r(T)dr(r) - 2r sa dSr(T). 



(B6) 
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The higher order corrections can be important at high 
temperature or energy. We also note that due to the pe- 
riodicity r(0) = r((3), Jq drrQpp = 0. Thus, at low energy 
the measure for the boson functional integral becomes 



where 



D[r,6] ~Y[r sa d6r{T)dQ(T) 



(B7) 



This is of a simple Cartesian product form, so that we 
don't have to deal with the interactions which might be 
generated by the nontrivial Jacobians of functional inte- 
gral measure. 

Now we turn to the details of the calculation of 
Eq. l(33Jl. Let us rewrite Eq. as follows: 

SS = -^2 Trln[(-G 1 )(l-G M)] 

rn 

= Trln[-Go *] " E Trln ( 1 ~ G ° M ) 

rn m 

~ -^TrH-Go^+^TrfGoM] 



\ ]T Tr [G MG M] . 



(B8) 



In the absence of the source field <j) m the sum over m 
would just give an overall factor N. The first term of 
the last line of Eq. (|B8|I is relevant for the mean-field ap- 
proximation. The second term ^2 m Trln[Go-M] contains 
a few parts. Among them gQ part is cancelled by the 
very definition of electron-phonon coupling, and the part 
linear in Sr is cancelled by the saddle point condition. 
The remaining nontrivial part is given by 

Trln [G M] f dr0 m G o 

[k}$ + K&] 8r(-iu), (B9) 



N 
~2 



where Krr and K$ are given by 



K<® = 2T^G ( l e)S (*e), (BIO) 

ie 

K$(iw) = 2T^Go(ie)[S (ie + zw)-E (ieX]Bll) 



Kg*(iu) = r 2 sa TY G (ie)Go(ie + iu) 



x (S (ic) + S (^e + ^w)) 2 , 
K&\iu) = -^Tj2Go(ie)G (ie + iuj), 

K re = T E G (ie)G (ie + icj) 



(B13) 



APPENDIX C: THE EXPLICIT FORMS OF THE 
POLARIZATION FUNCTIONS K (i \iuj) AT ZERO 
TEMPERATURE 

The polarization functions K(ioS) have been calculated 
in Appendix A of Ref. ^3 i n terms of the digamma func- 
tion at finite temperature. We will instead calculate the 
polarization functions at zero temperature in a closed 
form. The integrations are completely elementary, and 
only the results and the low frequency asymptotic behav- 
iors will be displayed. First we note that only the real 
part of Ki-r ■ (iu>) which is an even function of frequency 
contributes to the effective action. 



ReK$(iu) = — In 

7T 



n\w\ 



(M + A) s 



3" 



r 2 



(H + A) 2 + g 2 



A 2 + g 2 



(CI) 



(C2) 



The sum of the above two gives 
A U + 2A 



In 



7T \UJ\ 

And the rest are given by 
2A 



1+ ,o;l(H+2A) ( 



A 2 



tan 



ej + \ui\A 



(C4) 



The third term \ J2 m Tr ln[G MG M] can be re- 
expressed as follows: 

~ Trln [G MG M] = — ^ K&) 5r{iLo)8r{-iuo) 

5 E E *M r -( iQ + ^ + 

x(i6 + -^=Q + ^ m )(-iw) 
v iv 



4 E E 2i ^ } <M^) »-«,(ie + -j=Q + <M({m) 



A 



tt w (M + 2A) 



In 



H(2A+M) 



A 2 



The low (imaginary) frequency asymptotic behaviors are 
given by 



4A r A 2 A|w|e 



A 2 

K 



A 4 



v{x)u \ 2 Ao 



A 2 



(C6) 



(C7) 
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